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ABSTRACT 

We investigate the ability of single top quark production via qq' — > squark — > tb and 
qq' — > slepton — > tb at the LHC and Tevatron to probe the strength of R-parity violating 
couplings in the minimal supersymmetric model. We found that given the existing bounds on 
R-parity violating couplings, single top quark production may be greatly enhanced over that 
predicted by the standard model, and that both colliders can either discover R-parity violating 
SUSY or set strong constraints on the relevant R-parity violating couplings. We further found 
that the LHC is much more powerful than the Tevatron in probing the squark couplings, but 
the two colliders have comparable sensitivity for the slepton couplings. 
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1. Introduction 



The HERA data showed excess events in deep-inelastic positron-proton scattering at high- 
er 2 and high x, which are in apparent conflict with the Standard Model expectation [1]. The 
excess events have been interpreted as evidence of R-parity breaking supersymmetry [2] . Hence 
detailed examination of effects of R-parity breaking supersymmetry in other processes is in 
order. Some of the phenomenological implications of R-parity violating couplings at e + e~ 
colliders have been investigated in Ref . [3] . Constraints on the R-parity violating couplings have 
also been obtained from perturbative unitarity [4,5], n — n oscillation [5,6], z/ e -Majorana mass 
[7] , neutrino- less double j3 decay [8] , charged current universality [9] , e — /i — r universality [9] , 
— e scattering [9], atomic parity violation [9], deep- inelastic scattering [9], double nucleon 
decay [10], K decay [11,12], r decay [13], D decay [13], B decay [14-16] and Z decay at LEPI 
[17,18]. Another important effect of R-parity violating couplings is that they may enhance the 
flavor changing top quark decays to the observable level of the upgraded Tevatron and LHC 
[19]. 

As is shown in Refs. [20-24] , single top quark production is very interesting to study at the 
Tevatron and the LHC since, in contrast to the QCD process of it pair production, it can be 
used to probe the electroweak theory. Single top production processes have been used to study 
the new physics effects involving the third-family quarks in a model independent approach [25] 
and in specific models [26,27]. More recently, motivated by the evidence of R-parity breaking 
supersymmetry [28,29] from the anomalous events at HERA [1], single top quark production 
qq' — > tb at the Tevatron induced by baryon-number violating (BV) couplings A" (via the 
exchange of a squark in the t-channel) and by lepton-number violating (LV) couplings A' (via 
the exchange of a slepton in the s-channel) has been studied [30] in minimal supersymmetric 
model (MSSM). It was found [30] that the upgraded Tevatron can probe the relevant BV 
couplings efficiently, while the probe for the relevant LV couplings is very limited. 

In addition to the process qq' — ► tb mentioned above, which can be effectively studied at the 
Tevatron, the R-parity BV coupling can lead to the reaction qq' — > tb via an s-channel squark 
contribution which is suppressed at the Tevatron. This process is suitable for study at the LHC 
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and it probes a different set of BV couplings. In this paper we make a detailed study of the 
s-channel BV effect. For completeness we study the effect at both the LHC and the upgraded 
Tevatron. We also study the s-channel slepton contribution to qq' — > tb at the LHC which we 
compare to the result obtained for the upgraded Tevatron [30] . 

This paper is organized as follows. In Sec. 2 we present the Lagrangian for R-parity vi- 
olating couplings and squared matrix elements for the processes qq' — > squark — > tb and 
qq' — > slepton — > tb. In Sec. 3 we evaluate the signal for these processes and the SM back- 
ground, and give the probing potential of the LHC in comparison to the upgraded Tevatron. 



2. tb and tb production in R-parity violating MSSM 

2.1 Lagrangian of R-parity violating couplings 

The R-parity violating part of the superpotential of the MSSM is given by 

Wy = X ijk L i L j E c k + \' ljk L l Q 3 D c k + X» jk U?Dp e k + ^UH 2 . (1) 

Here Li(Qi) and Ei(Ui, Di) are the left-handed lepton (quark) doublet and right-handed lepton 
(quark) singlet chiral superfields, k are generation indices, and c denotes charge conjugation. 
Hi i2 are the chiral superfields representing the two Higgs doublets. The A and A' couplings 
violate lepton-number conservation, while the A" couplings violate baryon-number conservation. 
The coefficient A^ is antisymmetric in the first two indices and \'^ k is antisymmetric in the 
last two indices. In terms of the four- component Dirac notation, the Lagrangians for the A' and 
A" couplings that affect single top production at the Tevatron and the LHC are given by 



*144 + 44*4 + (4)*(4) c 4 



-e!44 - «£d*ei - (4m)X] + h.c, (2) 

£a» = -a^44(^)x + 4(4)X + ^K) c 4]+^-c.. (3) 

The terms proportional to A are not relevant to our present discussion and will not be con- 
sidered here. Note that while it is theoretically possible to have both BV and LV terms in 



the Lagrangian, the non-observation of proton decay imposes very stringent conditions on their 
simultaneous presence[31]. We, therefore, assume the existence of either LV couplings or BV 
couplings, and investigate their separate effects in single top quark production. 

2.2 qq' — > squark — > tb 

Production of tb via an s-channel diagram u' l d j — > d k R — > tb can be induced by the BV 
couplings A". The matrix element squared is given by 

2 

(Pi ' P2) [Ps ■ Pa ~ M t (s t ■ p 4 )\ , (4) 

where p\ and P2 denote the momenta of the incoming quarks u l and d? , P3 and P4 of the outgoing 
t and 6 quarks. The center-of-mass energy of the parton is given by s and s t denotes the spin 
of top quark which is given by 

s t = j^(\P3\,Etp 3 ), (5) 

where h — ±1 denotes the two helicity states, and ^3 is the unit three-vector in the momentum 
direction of top quark. 

Neglecting the contribution of third-family sea quark in the initial states, we obtain 

Y^\M x „(ud^tb)\ 2 

Y,\M x „{us^tb)\ 2 

Y^\M xl/ {cd^tb)\ 2 

Y^\M x „{cs^tb)\ 2 

In the R-parity conserving MSSM, the down-type squark d k R can decay into charginos and 
neutralinos via the processes d R ^u k + Xj (j = 1, 2) and d R ^d k + x°j (j = 1, 2, 3, 4), where 
Xj and x°j represent a chargino and neutralino, respectively [32]. Of course, it can also decay 
into gluino plus quark if kinematically allowed. In the R-parity violating MSSM, the down-type 
squark d k R can also decay into quark pairs d k — > d? + u l via the A" couplings. Since some of 
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the relevant A" couplings may be quite large, the width of a heavy down-type squark d k R can 
be large even if we do not consider the decay involving the gluino. We found that within the 
allowed parameter space (A", chargino and neutralino sector) its width Yg can be as large as 

2.3 qq' — > slepton —> tb 

Production of tb via an s-channel slepton u l d? — > e\ — > tb can be induced by the LV 
couplings A'. The matrix element squared is given by 

E! m«|* = iMg^l, • n) • * - M t(Sl • „„ , m 

where we assumed mass degeneracy for sleptons of different flavors. 

In the R-parity conserving MSSM, the charged sleptons e"x will decay into charginos and 
neutralinos via the processes e"x — > v e + Xj U = 1>2) and ex — > e + x° (j = 1,2,3,4) [32]. 
However, in the R-parity violating MSSM, the slepton can also decay into quark pairs via the 
A' couplings e l L — > u 3 L + d k R . Since the allowed A' couplings are quite small, the dominant decays 
are the chargino and neutralino modes. The partial widths are given by 

T(e L ^u e + ^) = ^q\U 3l f [Ml - Ml^)\ (11) 

r(e~x-e + x?) = ^|^i + ^-^ 2 | 2 (m|-M|o) 2 , (12) 

where sw = sin 9 w ,cw = cos9 w and the masses of the lepton and down- type quarks are 
neglected. The masses of charginos and neutralinos, and the matrix elements and N[- 
which respectively diagonalize the mass matrix of chargino and neutralino, depend on the 
SUSY parameters M 2 , M 1; and tan/3 [29]. Here, M 2 and M x are the masses of gauginos 
corresponding to SU(2) and U{\), respectively, /x is the coefficient of the HiH 2 mixing term in 
the superpotential, and tan/3 = v 2 /vi is the ratio of the vacuum expectation values of the two 
Higgs doublets. 



3. Numerical calculation and results 
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Due to the large QCD backgrounds at hadron colliders, it is very difficult, if not impossible, 
to search for the signal from the hadronic decays of the top quark. We therefore look for events 
with t — > W + b — > l + vb (I = e, //). (We take into account of the fact that the top quark is 
polarized in hadronic production.) Thus, the signature of this process is an energetic charged 
lepton, missing E?, and double 6-quark jets. We assumed silicon vertex tagging of the 6-quark 
jet with 50% effeciency and the probability of 0.4% for a light quark jet to be mis-identified as 
a fo-jet [33]. 

Although the present events have the unique signal of same sign b-quarks, since the tagging 
can not distinguish a b-quark jet from 6-quark jet, there are many potential SM backgrounds 
[33]: 

(1) the Drell-Yan like process qq' — > W* — > tb; 

(2) the quark-gluon process qg — > q'tb with a W-boson as an intermediate state in either 

the t-channel or the s-channel of a subdiagram; 

(3) processes involving a b-quark in the initial state, bq(q) — > tq'{q') and gb — > tW; 

(4) Wbb; 

(5) Wjj; 

(6) tt -> W~W + bb. 

Background process (2) contains an extra quark jet and can only mimic our signal if the quark 
misses detection by going into the beam pipe. This can only happen when the light quark jet 
has the pseudorapidity greater than about 3 or the transverse momentum less than about 10 
GeV. In our calculation of the W^-gluon fusion process as a background, we impose rj(q') > 3 
and Pt(q') < 10 GeV for the light-quark jet. The bq(q) — > tq'(q') background is greatly reduced 
by requiring double b-tagging. The process gb — > tW can only imitate our signal if the W 
decays into two jets, where one jet is missed by the detector and the other is mis- identified as 
a b quark, which should be negligible. Since we required two 6-jets to be present in the final 
state and assumed the probability for a light quark jet to be mis- identified as a 6- jet is 0.4%, 
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the potentially large background process (5) from Wjj is reduced to an insignificant level. Also 
we required the reconstructed top quark mass M (bW) to lie within the mass range 



\M(bW) - m t \ < 30 GeV, (13) 

which can also reduce the backgrounds Wbb and Wjj efficiently. Background process (6) can 
mimic our signal if both W's decay leptonically and one charged lepton is not detected, which 
we assumed to occur if rj(l) > 3 and pr(l) < 10 GeV. 

To make a realistic estimate we also need to consider the detector acceptance. To simulate 
the detector acceptance, we made a series of cuts on the transverse momentum (pr), the 
pseudo-rapidity (rj), and the separation in the azimuthal angle-pseudo rapidity plane ( AR = 
^/(A0) 2 + (Ar]) 2 ) between a jet and a lepton or between two jets. For the upgraded Tevatron, 
the cuts are chosen to be 

(14) 
(15) 

Vjj. U.V, • (16) 



For the LHC, the cuts are chosen to be 
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(17) 
(18) 
(19) 
(20) 
(21) 



To make the analyses more realistic, we simulate the detector effects by assuming a Gaussian 
smearing from the energy of the final state particles, given by: 

AE/E = 30%/ VE © 1%, for leptons , (22) 
= 80%/ VE © 5%, for hadrons , (23) 



where © indicates that the energy dependent and independent terms are added in quadrature 
and E is in GeV. 

We calculated the pp (for the Tevatron) and pp (for the LHC) cross sections for the sig- 
nal with the MRSA' structure functions [34]. We have also examined the effect of using the 
CTEQ3M [35] structure functions and found the difference between the two sets of structure 
functions to be small. We have explicitly calculated backgrounds (1) and (2), and for the oth- 
ers used the Wjj background analysis of Ref.[36]. The effect of the cuts is shown in Table 1. 
Also, in our numerical calculation, we assumed M t = 175 GeV, y/s = 2 TeV for the upgraded 
Tevatron and y/s = 14 TeV for the LHC. The integrated luminosities for both colliders are 
assumed to be 10 fb -1 . Assuming Poisson statistics, the number of signal events required for 
discovery of a signal at the 95% confidence level is approximately: 

73TT? £3 ' (24) 

where S (B) is the number of signal (background) events obtained by multiplying the signal 
(background) cross section by the luminosity ( 10 fb -1 ) and the tagging efficiency for two 6-jets 
(0.5 x 0.5). 

With all the above assumptions, we now present the results for both processes. 
3.1 The B- violating process of tb production 

For the BV process of tb production, qq' — > squark — > tb, we neglect ud — > s — > tb and 
us — > d — > tb since A" 12 < 10~ 6 [10]. For simplicity, we assume A^ and A^ do not coexist and 
hence evaluate cd — > s — > tb and cs — > d — > tb seperately. 

Assuming T s = Ms/5, we obtain Fig. 1 which shows the value of A^A^ versus strange- 
squark mass for cd — > s — > tb to be observable at 95% confidence level. The region above each 
curve is the corresponding observable region. The solid curve is for the LHC, the dotted curve 
is for the upgraded Tevatron and the dashed line is the perturbative unitarity bound [4,5]. 
Here we see that both the LHC and the upgraded Tevatron can efficiently probe the relevant 
couplings, and the LHC serves a more powerful probe than the upgraded Tevatron. 

As was discussed in the above section, the width of a down-type squark depends on many 



free parameters, which can vary in a large range. In order to show the sensitivity of the results 
to the width of strange-squark, we present in Fig. 2 the value of A^A^ versus the ratio Yg/Mg 
for cd — > s — > tb to be observable at 95% confidence level. Here we assume M s = 300 GeV. The 
region above each curve is the corresponding observable region. The solid curve is for the LHC 
and the dashed curve is for the upgraded Tevatron. We see from this figure that the value of 
A212A332 varies mildly as a function of Yg/Mg. Again it is shown in this figure that the LHC is 
more powerful than the upgraded Tevatron in probing the relevant couplings. 

The value of A^A^ versus down-squark mass for cs — > d — > tb to be observable at 95% 
confidence level is shown in Fig. 3. The region above each curve is the corresponding observable 
region. The solid curve is for the LHC, the dotted curve is for the upgraded Tevatron and the 
dashed line is the perturbative unitarity bound. The behaviour of this figure is similar to Fig.l. 
But for equal squark mass the value of A^A^ in Fig. 3 is higher than the value of A^A^ in 
Fig.l. This shows that the process cs — > d — > tb cannot be probed as efficiently as cd — > s — > tb 
because of the relative suppression of the strange quark structure function compared to the 
valence down quark. 

3.2 L-violating process of tb production 

For the LV process of tb production, qq' — > slepton — > tb, we only consider the dominant 
process ud slepton — > tb and thus provide the results for A / 111 A / 133 + A 211 A 23 3 + A^A^. The 
previous study [30] of this process at the upgraded Tevatron has shown that within the allowed 
range of the relevant coupling constants, this process is observable only when the slepton mass 
lies in a specific narrow range. Here we will determine if the LHC can do better than the 
upgraded Tevatron. 

As discussed in the above section, the allowed A' couplings which induce a charged slepton 
to decay into quark pairs are quite small and thus the dominant decays of a charged slepton 
are the chargino and neutralino modes. So we only consider the chargino and neutralino modes 
for simplicity. Then the width of the charged slepton only depends on the SUSY parameters 
M 2 , Mi, 11 and tan/3. In our calculation we use the GUT relation Mi = | a 2-M 2 ~ \Mi, and fix 
M 2 = —fx = 250 GeV and tan (3 = 2. We checked that in this case the chargino and neutralino 
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masses are above the present lower limits from LEP II[37]. 

Figure 4 shows the value of A' 111 A / 133 + A^A^ + A 311 A 3 33 versus the slepton mass for ud — > 
slepton — > tb to be observable at 95% confidence level. The region above each curve is the 
corresponding observable region. The solid curve is for the LHC, the dotted curve is for the 
upgraded Tevatron and the dashed line is the value obtained by considering the following bounds 
for squark mass of 100 GeV [7,11,17] 



\Kn\ 
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0.012, (i = 1,2,3), 


(25) 


1 V 1 

l A 133l 
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0.001, 


(26) 


l A 233l 


< 


0.16, 


(27) 


1 \' 1 

1 A 333 1 


< 


0.26. 


(28) 



Figure 4 shows that below the present upper limit for the couplings the LHC cannot do much 
better than the upgraded Tevatron in further probing the couplings. 

In summary, we have studied single top quark production via qq' — > squark — > tb and 
qq' — > slepton — > tb at the Tevatron and the LHC in the MSSM with R-parity violation. Our 
results show that from the measurement of single top production, the LHC can efficiently probe 
the relevant R-parity violating couplings. 
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Figure Captions 



Fig. 1 The value of A^A^ versus strange-squark mass for cd — > s —>■ tb to be observable 
at 95% confidence level. The region above each curve is the corresponding observable region. 
The solid curve is for the LHC, the dotted curve is for the upgraded Tevatron and the dashed 
line is the perturbative unitarity bound [4,5]. 

Fig. 2 The value of A^A^ versus the ratio Fg/Mg for cd — > s — > tb to be observable at 
95% confidence level. The region above each curve is the corresponding observable region. The 
solid curve is for the LHC and the dashed curve is for the upgraded Tevatron. 

Fig. 3 The value of A^A^ versus down-squark mass for cs — > d —>■ tb to be observable at 
95% confidence level. The region above each curve is the corresponding observable region. The 
solid curve is for the LHC, the dotted curve is for the upgraded Tevatron and the dashed line 
is the perturbative unitarity bound [4,5]. 

Fig. 4 The value of A / 111 A / 133 + A^A^ + A 311 A 3 33 versus the slepton mass for ud — > I — > tb 
to be observable at 95% confidence level. The region above each curve is the corresponding 
observable region. The solid curve is for the LHC, the dotted curve is for the upgraded Tevatron 
and the dashed line is the present bound, Eqs. (|25| - |28| ). 
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Table 1: 

Signal and background cross sections in units of fb after various cuts at the Tevatron and the 
LHC. The qq' — > q — > tb results have been calculated using the unitarity bound for the relevant 
couplings and assuming Mg = 500 GeV and = Mg/5. The qq' — > / — ► tb results have been 
calculated using the present upper bound Eqs . (|25|) - (|28|) for the relevant couplings and slepton 
mass of 300 GeV. The slepton width is calculated by assuming M 2 = — \i = 250 GeV and 
tan/3 = 2. The charge conjugate channels are included. 



Tevatron 


basic cuts 


basic+m(bW) cuts 


basic+m( b Wj + bb-tag 


cd — > s — > tb 


1545 


1436 


359 


cs — > d — > tb 


186 


174 


44 


ud — > / — >• tb 


75 


73 


18 


qq' — > tb 


78 


75 


19 


gq -> q'tb 


4 


3.4 


0.85 


qb — > q't 


236 


224 


0.45 


Wbb 


264 


122 


30 


Wjj 


62900 


45000 


0.7 


tt 


16 


7 


1.8 


LHC 


basic cuts 


basic+m(bW) cuts 


basic+m(bW)+bb-tag 


cd — > s — > tb 


335600 


304800 


76200 


cs ^ d tb 


125200 


113600 


28400 


ud — > / — > t6 


482 


480 


120 


gg' — > t6 


573 


547 


137 


5-g -> g'tfe 


1104 


810 


203 


g6 — >■ q't 


14150 


13440 


27 


Wbb 


756 


338 


85 


Wjj 


623600 


379000 


6 


tt 


1840 


644 


161 
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II II II 

^11^33+^811^233+^311^333 




